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Electromagnetic interactions are incorporated into Heavy Hadron Chiral Perturbation 
Theory. Short and long distance magnetic moment contributions to the chiral Lagrangian 
are identified, and Ml radiative decays of heavy vector mesons and sextet baryons are 
studied. Using recent CLEO D* branching fraction ratio data, we fit the meson coupling 
to the axial vector Goldstone current and find gf = 0.34 ±0.48 for rric = 1700 Me V. Finally, 
we obtain model independent predictions for total and partial widths of charm and bottom 
vector mesons. 
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A synthesis of Chiral Perturbation Theory and the Heavy Quark Effective Theory 
(HQET) has recently been developed [jl|-|5|. This hybrid effective theory describes low 
energy strong interactions between light Goldstone bosons and hadrons containing a single 
heavy quark. Weak b c transitions among heavy meson or baryon states can also be 
incorporated into this framework. In this letter, we extend the theory's formalism to 
include electromagnetism and then study the radiative decays of heavy vector mesons and 
sextet baryons. 

To begin, we briefly review the basic elements of Heavy Hadron Chiral Perturbation 
Theory (HHCPT). The Goldstone bosons resulting from the chiral symmetry breakdown 
SU{3)i^ X SU{3)ji 5't/(3)^_|_H appear in the pion octet 



8 1 

a=l 



(1) 



and are associated with the pion decay constant / ^ 93 MeV. These fields are arranged 
into the exponentiated matrix functions E = e^'^'f and ^ = "v^" = e'^'f that transform 
under the chiral symmetry group as 

(2) 

Here L and R represent global elements of SU{3)l and SU{3)in while U acts like a lo- 
cal SU{3)i^-^-ji transformation. Chiral invariant terms that describe Goldstone boson self 
interactions are constructed from the fields in (0) and their derivatives. 

Hadrons containing a heavy quark emit and absorb light Goldstone bosons with no 
appreciable change in their four velocities. They are consequently described by velocity 
dependent fields. In the meson sector, we introduce the operators Pi{v) and Pi^iy) that 
annihilate pseudoscalar and vector mesons with quark content Qq. When the suppressed 
heavy quark label carried by these fields corresponds to charm, their individual components 
are given by 

{P^,P2,P^) = {D^,D+,Dt) 

(3) 

{p:,p;,p;) = {D*\D*+,D*t). 



^ This introductory discussion closely follows that presented in refs. [^J5[ to which we refer 
interested readers for further details. 
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In the infinite quark mass limit, it is useful to combine the degenerate meson spin states 
into the 4x4 matrix field 

H^iv) = ^ [-P.ivh' + P:,{v)r] (4a) 

and its conjugate 

ifiv) = joH^o = [p^\vh' + p;^\v)r] m 

H carries a heavy quark spinor index and a separate light antiquark spinor index and 
transforms as an antitriplet under SU{?>)lj^r and doublet under SU{2)y. 

Baryons with quark content Qqq enter into the theory in two types depending upon the 
angular momentum of their light degrees of freedom. In the first case, the light spectators 
are arranged in a symmetric spin-1 configuration that couples with the heavy spin-i quark 
to form = and = f ^ states. When the heavy partner is taken to be charm, 
the spin- 1 states are destroyed by the Dirac operators appearing in the symmetric sextet 
representation 



1=1 



iv+ to /]~^o' 

2^c Y 2"c 



(5) 



\ Y 2"c \l 2^c -^'c / 

Their spin-| counterparts are annihilated by the corresponding Rarita-Schwinger field S*. 
We again combine the Dirac and Rarita-Schwinger operators into the "super" fields 0] 

_ (6) 

Then 5"^ transforms as a sextet under SU{3)i^+uj doublet under SU{2)y, and is an axial 
vector. 

The spectators in the second case are bound together into an antisymmetric spin-0 
state. The resulting = baryons are assigned to the field Ti{v), which is an SU (3)^+^ 
antitriplet and SU (2)„ doublet. When Q = c, the components of Ti are the singly charmed 
baryons 

m,T2,T3) = (S«,-S+,A+). (7) 



These antitriplet baryons can alternatively be arranged into the antisymmetric matrix 
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-^ -A+ 



2"c 



2 c 



1^0 
2"c 



2"c * 



1^0 
2"c 



(8) 



where {T-^^) = eijk/V^- 

We can now construct the zeroth order effective chiral Lagrangian that describes 
the low energy interactions between light Goldstone bosons and heavy hadrons in the 
infinite heavy quark mass limit. The leading order terms must be hermitian, Lorentz 
invariant, light fiavor and heavy quark spin symmetric, and parity even. We can also 
readily incorporate electromagnetism into the hybrid chiral theory by gauging a U{1)em 
subgroup of the global S't/(3)^ x 5'C/(3)h symmetry group. Only long wavelength photons 
with energies less than the chiral symmetry breaking scale explicitly remain in the low 
energy theory while short wavelength modes are integrated out. In d = 4 — e dimensions, 
the effective Lagrangian looks like i 



£(0) 
£(0) 



-^F^^'F^, + ^i^TV(P^Etp^E) (9a) 

E {-iTr{H'\ ■ VH[} - i%v ■ VS'J + {Ms - Mt)%S'^ + iT\ -VT, 
Q=c,b 



+ giTr{Hl{4))j'H")+ig2e^..xS^,v-'{A'^)]{S 



X\jk 



+ 93 



(96) 



The Goldstone bosons explicitly couple to the matter fields through the axial vector com- 
bination 



(10a) 



They also communicate via the vector field 



2 



(106) 



^ Meson contributions are written in terms of the dimension- 1 field H' = \/ MhH so that all 
heavy mass dependence is removed from the leading order Lagrangian. 
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that appears inside the heavy hadron covariant derivatives □ 

V^S'J = d^S^ + {V^y^^Sl^ + {y^'MSf - i^'/'^eA^ [QqS'J + Q^^Sl^ + 2^5^^] (11) 



The remaining Goldstone covariant derivatives in and (|T^) are given by 



(12) 



where 

2=1 22 I = r I (13) 




3 

denotes the light quark electric charge matrix. 

Spin symmetry violating contributions to the chiral Lagrangian enter at ©(l/mg). 
Among these are heavy quark magnetic moment terms which mediate Ml radiative tran- 
sitions. As we will see, these terms are completely fixed by heavy quark number conser- 
vation. This simple but crucial observation allows one to use experimental meson decay 
information to determine the parameter gi. 

Recall that the photon gauge field couples to the conserved current that counts heavy 
quark number in the underlying QCD theory as well as in the low energy HQET and 
HHCPT. The original QCD current appears in its well-known Gordon decomposed form 
as 

Running down in energy to the heavy quark thresholds and invoking the velocity super- 
selection rule to set p^'^ = niQV + k^'\ one can match this tree level current onto the 
corresponding HQET expression 

^ J^ic^) + ^0,- t,) + ^CT,.(^^ + hi''^ (15) 



^ We distinguish the photon field A'^ from the axial vector Goldstone current A*^ by writing 
the former in calligraphy type. Similarly, we let Q represent the electric charge operator, which 
is different from the heavy quark symbol Q. 
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The gauge field couples to this current in the Lagrangian 



(16) 



Q=c,b 



where the 0{l/mQ) Oi operators are constructed from either symmetric or antisymmetric 
combinations of two HQET covariant derivatives = — ifi'^/'^ gG^Ta — i^^^'^eA^Q: 



Oi 



2m. 



02 = ^/^i^V^.T^/ii'^^G'r 



4m 



Q 



^^/'eQQ -(g) (Q) 



(17a) 
(176) 
(17c) 



The ai coefficients of these dimension-five operators are thus fixed by (|T5|) and equal unity 
at lowest order J^. 

Running down further in energy from the heavy quark thresholds to the chiral sym- 
metry breaking scale A^, we match the electromagnetic pieces of operators Oi and O3 onto 
the following short distance contributions to the HHCPT Lagrangian: 



^ishort) =Y,{- 



Q=c,b 



2m 



Amr 



2mQ 
1 

2mc 



4mr 



^'/^QQe(mg)-^ 



Amr 



1- 



(18) 



These 0{l/mQ) terms describe the interaction of photons with the heavy quark constituent 
inside a if', S or T hadron. Consequently, the Lorentz and fiavor indices for the hadrons' 
light degrees of freedom are trivially contracted. Heavy quark number conservation deter- 
mines the ratio of the operator coefficients in (p!8|) to the kinetic terms in the zeroth order 
Lagrangian (|9^). The conserved HQET current therefore matches onto 



JHHCPT 



1 1 



TrH 



Vn + 



2mQ 



2m, 



H' 



+ T 



2m. 



2m. 



2mr 



(19) 



in the low energy chiral theory. 

Photons also couple to the light brown muck inside heavy hadrons leaving the spins 
of their heavy quark constituents unaltered. Such long distance interactions generate 
additional electromagnetic contributions to the effective Lagrangian at the scale. We 
focus upon just the induced magnetic moment terms: 



(long) 



A. 



(20) 



A few points about these long distance operators should be noted. Firstly, the suppressed 
heavy quark spinor indices in these spin symmetry preserving terms are simply contracted. 
Their light Lorentz and flavor indices on the other hand are nontrivially arranged. Sec- 
ondly, the coefficients Ch, Cs and Cst are a priori unknown. But naive dimensional analysis 
suggests that they are of order one fl^. Finally, there is no long distance magnetic mo- 
ment interaction for just the antitriplet baryon since the photon field cannot couple to its 
spinless light degree of freedom. 

Having identified the short and long distance magnetic moment terms in the low energy 
chiral theory, we can now study Ml radiative transitions between meson and baryon states. 
Since the hyperfine splitting between charmed pseudoscalar and vector meson partners is 
only slightly greater than a pion mass, the electromagnetic decay D* D'j competes with 
the strong process D* Dn. The greater phase space for the electromagnetic transition 
offsets its inherently smaller amplitude. Bottom vector mesons must radiatively decay 
because pion emission is kinematically forbidden. So these Ml meson processes are of 
genuine phenomenological interest. Similar considerations apply to the baryon transitions. 

The vector meson and sextet baryon radiative decay rates are readily determined from 



the magnetic moment terms in (18) and (|2 



FiP* -> Pa) = 



(21a) 



- - 



M, 



3 VMp 

r(^*' ^ ^^7) = 

1 /M, w M|.-M| x3 
ISVMsJv Ms* ) 

r(^(*)^^T,7) = 

1 / Mr ^/^sV) -^t\3 



(216) 



Qc 



-aEM[mQ, 



)V2 + 2^IV(T^^,/QT^^,>,^(A^)i 



L/2 



f21c) 



6 VM 



c 

LA 



X 



One can clearly identify the short and long distance contributions to these partial widths 
from their electric charges and associated inverse mass scales. The corresponding strong 
interaction decay rates are derived from the Goldstone axial vector couplings in the leading 
order Lagrangian 
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487r/2 VMp* 



9i 



[Ml, - (Mp + m^)^][M|. - (Mp - m^y 
Mt 



[M|, - {Ms + m^)^] [Ml - {Ms - 
Mi 



3/2 



J* I 



3/2 



(22a) 



(226) 



(22c) 



9l 



[M2(., -(M, + m^)2][M2., 



(M, 



Tot, 



3/2 



|Tr(Tf V"T(^,,)|2. 



247r/2 VM<,(*) 

None of the heavy hadron electromagnetic and strong partial widths have been directly 



measured. However, values for D* branching fraction ratios are known [|TT 



P^ 



D+ _ 

^7 p^^* + 



r(D*+ ^ D+7) 



D+ttO) 

Taken in conjunction with the isospin relation 



0.572 ±0.057 ±0.081 
= 0.035 ±0.047 ±0.052. 



(23) 



Rt 



these data yield the following branching fractions: S 



2.21 ±0.07, 



(24) 



]J*+ _ 




68.1 ±1.0 ±1.3% 


(25a) 


]J*+ _ 




30.8 ±0.4 ±0.8% 


(256) 




D+7 


1.1 ± 1.4 ± 1.6% 


(25c) 



63.6 ±2.3 ±3.3% (25d) 
D*o^DO^ 36.4 ±2.3 ±3.3%. (25e) 

^ These very recent CLEO values differ significantly from Particle Data Group world averages 
p. 
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Using the branching fraction ratios for the two independent D* charge modes in (p3[), 
we can deduce the parameters Ch/-^x 9i that enter into the heavy meson electromag- 
netic and strong decay rates respectively. To extract these unknown couplings from the D* 
data and to predict the B* widths, we must specify numerical values for the charm and bot- 
tom mass parameters rric and mf, . Since these quark masses are sources of large theoretical 
uncertainty, we perform the fit twice. First we assume (mcTOb) = (1500 MeV, 4500 MeV), 
and then we take [rric, mi,) = (1700 MeV, 5000 MeV). Reasonable estimates for the heavy 
quark masses are covered by the range between these two sets of input values. 

From the charm vector meson ratios, we find two equations for the two unknowns: 

2 



9i 



1 



2mr 



-a, 



1/2 _ 

A. 



a. 



r(A, 



,1/2 



(M2,+ -(M^++m^o)2)(M 



(M, 



m 



3/2 



(M 



9i 



2 

D* 
-I 2 



R. 



+ 



Arrir 



-a, 



I Nl/2 I '-H 



(26) 
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(M2 ,0 - (Moo + m,o)2) [Ml,, - (M_ 



DO 



3/2 



5127r/2 



[Ml 



Following the suggestion of naive dimensional analysis, we choose the roots of these 
quadratic equations that yield values for gf of order unity. The results of the parame- 
ter fit are then listed as functions of the charm quark mass in Table 1: 



Coupling 


rUc = 1500 MeV 


rric = 1700 MeV 


9l 


(-0.68±0.50)/(1000 MeV) 
0.43 ±0.61 


(-0.60 ± 0.44)/(1000 MeV) 
0.34 ±0.48 



Table 1 



These results only weakly depend upon A^^ through the logarithmic running of the fine 
structure constant. Therefore, a very precise numerical value for the chiral symmetry 
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breaking scale need not be specified. However, if one reasonably assumes ^ 1000 MeV, 
then the value for turns out to be of order one and is consistent with our earlier ex- 
pectations. We also note for comparison that the nonrelativistic quark model estimate for 
the squared Goldstone axial vector parameter is 0.7 <gf < 1.0 P,|l3l. The HHCPT central 
value for this coupling is therefore of the same order of magnitude but smaller than the 
quark model number. The large error bars on gf reflect the 200% uncertainty in the mea- 
surement ( |25cj) of the D*~^ — > D'^'j branching fraction. Improvements in the experimental 
value will yield more precise estimates for this basic chiral Lagrangian parameter. 

Having found c^ij and (7^, we can now obtain model independent predictions for the 
total and partial widths of all D* and B* vector mesons. Our predictions are summarized 
in Table 2: 
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Width ( MeV) 


nic = 1500 MeV 
nib = 4500 MeV 


nic = 1700 MeV 
rrib = 5000 MeV 


T{D*+) 

T{D*+ L»°7r+) 
T{D*+ D+^) 


(12.44 ±12.27) X 10-2 
(3.56 ±5.06) X 10-2 

(7.83 ± 11.13) X 10-2 
(1.06 ± 1.05) X 10-2 


(9.70 ±9.56) X 10-2 
(2.77 ±3.94) X 10-2 
(6.10 ±8.68) X 10-2 
(0.83 ±0.81) X 10-2 




(6.49 ± 7.94) X 10-2 
(5.36 ±7.63) X 10-2 
(1.13 ±2.20) X 10-2 


(5.06 ±6.19) X 10-2 
(4.18 ±5.94) X 10-2 
(0.88 ± 1.71) X 10-2 




(8.46 ± 11.94) X 10-^ 


(6.60 ±9.31) X 10-^ 




(1.63 ±2.61) X 10-^ 


(1.27 ±2.03) X 10-^ 



Table 2 



Current upper bounds on D* widths are about an order of magnitude greater than the 
central values quoted here, while no B* decay information is yet available. Comparison of 
these theoretical results with experimental data must therefore be left for the future. 

To conclude, we comment upon several possible extensions of this work. In the meson 
sector, a number of refinements of our leading order analysis should be pursued. Perturba- 
tive QCD corrections, subleading 0{l/mQ) and SU{3)l+r breaking effects, and calculable 
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nonanalytic terms from Goldstone boson loop diagrams may all be systematically incorpo- 
rated into the HHCPT framework to yield improved values for the meson parameters and 
decay rates. D* and B* decays can also be worked out and studied in a straightforward 
fashion. For the sextet baryons, the present absence of branching ratio data precludes 
our determining the baryon couplings (cs/Aj^,(7|) and {cst/-^x^ di) ^^^^ widths 
of the spin- 1 states in precisely the same manner as their meson analogues. Nonetheless, 
such baryon data will eventually become available in the future. So the enhancements 
mentioned above for the mesons ought to be carried out for the baryons as well. Finally, 
the scope of HHCPT can be broadened to include higher resonances such as the Di and 
D2 states Electromagnetic interactions for these meson and baryon excitations may 



be incorporated into the theory along the same lines as those for the heavy hadron H' , S 
and T ground states. 
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